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ABSTRACT 

Concern about systematic uncertainty in the 4 He abundance as well as the chemical evo- 
lution of 3 He leads us to re-examine this important limit. It is shown that with conservative 
assumptions no more than the equivalent of 4 massless neutrino species are allowed. Even 
with the most extreme estimates of the astrophysical uncertainties a meaningful limit still ex- 
ists, less than 5 massless neutrino species, and illustrates the robustness of this argument. A 
definitive measurement of the deuterium abundance in high-redshift hydrogen clouds should 
soon sharpen the limit. 



Introduction. Big-bang nucleosynthesis is one of the experimental pillars of the standard 
cosmology 0> @|. It also probes particle physics. Among other things, it has been used to set 
a stringent limit to the energy density contributed by light (mass <C MeV) particle species, 
usually quantified as the equivalent number of massless neutrino species (= N u ) ||. This 
limit indicated that the number of neutrino species was small before accelerator experiments 
were able to experimentally establish this directly, and further, has served to constraint 
theories proposed to unify the forces and particles of Nature. Since the cosmological bound 
also constrains light species that do not couple to the Z°, it is an important complement to 
the LEP measurement, N v — 2.991 ± 0.016. 

The physics underlying the neutrino limit is simple: the big-bang production of 4 He 
increases with both baryon density (quantified by the present baryon to photon ratio rj) 
and the number of massless neutrino species. Thus, an upper limit to the primeval 4 He 
abundance (= Yp) and a lower limit to the baryon density lead to an upper limit to N u . 
The lower limit to the baryon density is based upon the big-bang production of deuterium, 
which rises rapidly with decreasing baryon density |4| . Because D is easily destroyed in stars, 
first being burnt to 3 He, the limit actually hinges upon D + 3 He, bringing in the chemical 
evolution of 3 He. This has been the standard approach for setting a limit to N v . 

Over the past five years limits to N v ranging from 3.04 to around 5 have been quoted 
||. The disparity arises because the irreducible uncertainties are systematic, rather than 
statistical. In particular, the pressing issues involve the primeval abundance of 4 He and the 
chemical evolution of 3 He (astronomers use the term "chemical evolution" to refer to nuclear 
processing). The purpose of this paper is to clarify the current situation, to use a new 
technique to obtain a limit to N u that is independent of the chemical evolution of 3 He, and 
to show how measurements of the primeval D abundance in high-redshift hydrogen clouds 
should soon sharpen the neutrino limit. In so doing, we will emphasize the robustness of the 
cosmological limit to N u . 

The light- element abundances. To orient the reader we begin with a brief overview. The 
big-bang production of D, 3 He, 4 He and 7 Li is summarized in Fig. 1. The predicted and 
measured abundances are consistent — within the uncertainties — for 77 ~ (2 — 7) x 10~ 10 
|T|, The status of the four light elements is as follows. 

1. The lithium abundance is measured in the atmospheres of the oldest stars in our galaxy 
(pop II halo stars) and is relatively well determined, ( 7 Li/H) pop n = (1.4 ±0.3) x 10~ 10 
[P. [I. Because 7 Li burns at a low temperature, the most important concern is that 7 Li 
may have been depleted. Stellar models indicate that a depletion of up to a factor of 
about two is consistent with the observations of 7 Li and other light elements in these 
old stars |§. 

2. The 4 He abundance increases with time, as 4 He is made by stars, and so any mea- 
surement provides an upper bound to the primeval component. The primeval abun- 
dance is extrapolated from measurements of metal-poor, extragalactic HII regions (hot 
clouds of ionized H and He gas). A compilation and analysis of the extant data gave 



1 



Yp = 0.232 ± 0.003 (stat) ± 0.005 (sys) 0. Others have estimated the systematic un- 



certainty to be larger, by a factor of two or even more [ IP] ; it arises from many sources, 
from uncertainties in theoretical emissivities to assumptions and approximations made 
in modeling the HII regions. A very recent study of 27 new extragalactic HII regions, 



using a new set of emissivities, finds Yp = 0.243 ± 0.003 |TT 



3. Deuterium is destroyed in the contemporary Universe, being burnt to 3 He, which links 
its post-big-bang history to that of 3 He. The deuterium abundance has been measured 



in the local interstellar medium (ISM) and in the pre-solar nebula |T5]. While neither 
value is expected to reflect the primeval abundance, both place a lower limit to the 
big-bang production and thereby a firm upper limit to the baryon density. This two- 
decade old limit, flp ^ 0.06 (70 km sec -1 Mpc -1 /ifo) 2 , precludes baryons from closing 
the Universe [[13[], and is not being questioned. 

The primeval deuterium abundance can be determined by measuring the deuterium 
Ly-a absorption feature in high-redshift (z > 2) hydrogen clouds which are backlit 
by distant QSOs ]14|, H|. There are now three claimed detections [16| and four ten- 
tative detections flTfl; the reported values are from (D/H)p ~ 2 x 10~ 5 to 2 x 10~ 4 , 



spanning the range anticipated on the basis of the big-bang production of the other 
light-elements. A definitive determination will peg the baryon-to-photon ratio to an 
accuracy of around 10% because of the rapid variation of D production with r\. 

4. Since all measurements of 3 He involve samples that have to some extent been processed 
through stars, chemical evolution is the central issue. The chemical evolution of 3 He is 
complicated and not well understood. Stars burn D to 3 He before reaching the main 
sequence; this is not controversial. According to conventional stellar models, massive 
stars (M > 8M ) reduce the 3 He in material they return to the ISM, by a factor 
of around five or so. Since high-mass stars also produce elements heavier than 4 He 
("metals"), there is a limit to the amount of processing through high- mass stars - and 
thus to 3 He destruction. Low-mass stars (M < 2M Q ) preserve and add significant 3 He 
to the material they return to the ISM ||18|| . This leads to the conclusion that the 
sum of D + 3 He should have decreased by at most a modest factor since primordial 
nucleosynthesis, which has been used to obtain the lower bound to r\ needed to constrain 
the number of massless neutrino species jp. 

There is little empirical evidence for the conventional view concerning the chemical 
evolution of 3 He; moreover, it has recently been challenged both theoretically and 
observationally. Slow-mixing mechanisms that would transport 3 He deep enough to 
be burnt have been suggested to explain oxygen and carbon isotopic ratios measured 
in some stars and meteorites Hl9fl . In addition, even if the conventional wisdom is 



correct in the mean, galactic chemical abundances are heterogeneous and the D + 3 He 
abundance has only been measured locally |2(| . 



On the observational side, a recent measurement of the He abundance in the ISM 
indicates that the sum of D + 3 He has not increased in the past 4.5 Gyr [^1], p2f ; 
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conventional models predict an increase since the recent evolution of the ISM should 
be dominated by low-mass stars. Finally, if the highest values of the D abundance 
measured in high-redshift hydrogen clouds are correct, then the present D + 3 He 
abundance is too small by a factor of almost ten to accommodate conventional models. 

A new approach. The two major obstacles to setting a reliable limit to N v are the 
chemical evolution of 3 He - which affects the lower bound to 77 based upon D + 3 He - and 
the systematic uncertainty in the primeval 4 He - which affects the upper bound to big-bang 
4 He production. To clarify and to minimize the dependence upon these we have computed 
limits to N u based upon a likelihood function for N = N v — Ay/0.016 which uses the 7 Li 
and 4 He abundances^ and takes deuterium into account only by means of a Bayesian prior. 
This approach exploits the fact that 7 Li increases with N u rapidly enough to provide a 
meaningful lower limit to rj as well as the fact that N v and a systematic shift (AY) in the 
assumed primeval 4 He abundance are equivalent (see Ref. ||). 

In computing C(N) we take Yp = 0.242 + AY ± 0.003 and consider two possibilities 
for 7 Li, ( 7 Li/H) P = (1.4 ± 0.3) x 10~ 10 (no depletion) and (3 ± 0.6) x 10~ 10 (factor of two 
depletion), using whichever gives the less stringent limit to N v . Regarding the primeval 4 He 
abundance; as a "standard case" we take the central value Yp = 0.242, 0.01 larger than 
the central value recommended in Ref. || - essentially equal to the value derived from the 



newest data set [|Tl|] - and roughly equal to the upper limit to Y P used previously in deriving 
a limit to N„. Through AY we allow for possible systematic error in the 4 He abundance and 
explore the sensitivity of the neutrino limit to it. 

We consider four priors for deuterium: (1) (D/H)p < 1.0, corresponding to essentially 
no prior information; (2) (D/H)p < 2 x 10 -4 , corresponding to the previous bound based 
upon D + 3 He H|; (3) a distribution for (D/H) P based upon an extreme chemical evolution 
model [p^l , which allows for destruction of 3 He in low-mass stars and extra destruction in 
high-mass stars by assuming that metals they produce were ejected into the intergalactic 
medium, avoiding the Galactic metallicity constraint; and (4) (D/H)p = (2.5 ±0.75) x 10~ 5 , 
an example of the accuracy to which the primeval deuterium might be determined. Our four 
likelihood functions are shown in Fig. 2. 

Discussion. Credibility intervals for N u are obtained from £(N) by choosing a value for 
AF and imposing a prior on N„, either N v > 3 or N v > 2. As Table 1 illustrates, better 
knowledge of the deuterium abundance and/or a tighter limit to the primeval 4 He abundance 
improve the limit to N v . 

Because our likelihood is a function of iV = N v — Ay/0.016, in the absence of a prior for 
N v the limit to N v would scale precisely as the limit to N plus Ay/0.016; e.g., increasing Yp 
by 0.01 would raise the limit by about 0.6 neutrino species. This scaling works reasonably 
well for the N u >2 prior. 

Regarding the dependence of the neutrino limit upon deuterium: 



1 The authors of Ref. |23 have advocated the use of 7 Li and 4 He alone, albeit in a different context. 
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D/H ( 7 Li) 


Y P 


N v >2 


N v >3 




232 ± 003 


2 0-2 6 2 8-3 6 


3.0-3.7 




237 ± 003 

V / • '.J 1 — 1 — V / • V / V 7 '.J 


2 0-2 9 3 2-3 9 


3.0-4.0 


< 1 (high) 


0.242 ± 0.003 


2.3-3.3, 3.4-4.3 


3.0-3.3, 3.4-4.3 




0.247 ± 0.003 


2.6-3.6, 3.8-4.6 


3.0-3.6, 3.7-4.6 






9 q Q q 4 1 4 q 


3 D 3 q 4 1 4 q 




0.232 ± 0.003 


2.0-3.1 


3.0-3.3 




0.237 ±0.003 


2.3-3.4 


3.0-3.5 


< 2 x 10~ 4 (low) 


0.242 ± 0.003 


2.6-3.7 


3.0-3.7 




0.247 ±0.003 


2.9-4.1 


3.0-4.0 




0.252 ±0.003 


3.2-4.4 


3.2-4.4 




0.232 ± 0.003 


2.0-2.6 


3.0-3.5 




0.237 ±0.003 


2.1-2.9 


3.0-3.6 


extreme chemical 


0.242 ± 0.003 


2.4-3.2 


3.0-3.4 


evolution for 3 He (high) 


0.247 ±0.003 


2.6-3.7 


3.0-3.6 




0.252 ±0.003 


2.9-4.0 


3.0-3.8 




0.232 ± 0.003 


2.0-2.4 


3.0-3.1 




0.237 ±0.003 


2.0-2.7 


3.0-3.2 


(2.5 ± 0.75) x 10~ 5 (high) 


0.242 ±0.003 


2.3-3.0 


3.0-3.2 




0.247 ±0.003 


2.6-3.3 


3.0-3.4 




0.252 ±0.003 


2.9-3.6 


3.0-3.6 



Table 1: The 95% credible intervals for N u based upon the different deuterium and N„ 
priors and central values of the primeval 4 He abundance. High (low) indicates high (low) 
7 Li abundance was used to derive the limit (see text). 
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1. Using the original D + 3 He bound, (D/H)p < 2 x 10 4 Q, the 95% credible region 
extends to N = 3.3. For Y P = 0.242, N v < 3.7(3.7) and for Y P = 0.252, JV V < 4.4(4.4). 
(The numbers in parentheses are for the prior N v > 2.) 

2. Using the probability distribution for deuterium based an extreme model for the chem- 
ical evolution of 3 He, the 95% credible region only extends to N = 3.2. For Yp = 0.242, 
N v < 3.4(3.2) and for Y P = 0.252, N v < 3.8(4.0). The use of a probability distribution 
for (D/H)p actually improves the bound slightly, in spite of the extreme model for 
chemical evolution assumed. 

3. If (D/H)p is determined to be (2.5 ±0.75) x 10~ 5 from high-redshift hydrogen clouds, 
the 95% credible limit region extends to N = 3.0. For Y P = 0.242, N v < 3.2(3.0) and 
for Y P = 0.252, N v < 3.6(3.6). 

4. With no information about deuterium, (D/H)p < 1.0, the likelihood function for N is 
double-peaked. The two peaks correspond to two ways of accommodating both 4 He 
and 7 Li: low N and high rj and high N and low rj. The upper boundary of the 95% 
credible region is N = 4.3. For a central value Yp = 0.242, N„ < 4.3(4.3) and for 
Yp = 0.252 (four times the systematic error estimated in Ref. ||), N v < 4.9(4.9). 

Several points should be noted. (1) Because 7 Li depletion in pop II halo stars cannot be 
proved nor disproved at the moment, our limits are based upon the 7 Li abundance that gives 
the less stringent limit. Once the issue of depletion is resolved, the limits may improve by a 
modest amount. Our limits are relatively insensitive to even larger Li depletion. For example 
going from a factor of two depletion to a factor of three depletion changes the neutrino limit 
by at most 0.2. (2) If the primeval deuterium abundance is determined to be at the high 
end of the current range, say (D/H)p = (1.0 ± 0.3) x 10~ 4 , the limits are essentially as in 
case (1) above. This is because such a measurement would pin 77 at very nearly the same 
value as the D + 3 He bound. However, should (D/H)p ~ (2 ±0.6) x 10 -4 , the limit worsens 
slightly (to 4.3 for a central Y P = 0.242), because 77 is pinned to a lower value. (3) The limits 
quoted here are not directly comparable to other work. The most similar approach is that 
of Fields et al. p|; however, they determine a lower limit to r) based upon 4 He and 7 Li alone 
(assuming N u — 3), and then use it to set an upper bound to N v (as opposed to determining 
a likelihood function for N„). 

The Table illustrates the robustness of the cosmological limit to the number of light 
particle species. With conservative assumptions about the astrophysical uncertainties - 
the extreme model for 3 He chemical evolution and Yp = 0.252 - less than four massless 
neutrino species are allowed. A comparable limit follows with Yp = 0.247 ± 0.003 and 
(D/H)p < 2 x 10~ 4 or no bound to deuterium and Y P = 0.242 ± 0.003. Even with the most 
extreme assumptions - no limit to primeval deuterium and Yp = 0.252 ±0.003 - a meaningful 
limit still exists, N v < 4.9. Thus, until the current astrophysical uncertainties are clarified 
we would argue that N v < 4 is a conservative limit, and N v < 5 is a very cautious limit. 

Astrophysical and cosmological limits will always have irreducible systematic uncertainty 
that cannot be quantified with a standard error. In the example at hand, the correctness 
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of the standard model of nucleosynthesis must be assumed - though there is no reason to 
doubt or a compelling alternative - and abundances in the contemporary Universe must be 
extrapolated to their primeval values. The primary concerns that motivated our reexamina- 
tion of the big-bang neutrino limit are the systematic uncertainties in the primeval D and 
4 He abundances, with the former dominated at this time by the chemical evolution of 3 He. 
We have shown that these concerns do not preclude setting a meaningful and robust limit to 
N v . With conservative assumptions, consideration of 7 Li and 4 He alone provides the limit, 
N v < 4. Even with extreme assumptions, the limit N v < 5 follows. The uncertainty due to 
the chemical evolution of 3 He will vanish when the primeval abundance of D is determined 
unambiguously in high-redshift hydrogen clouds, which, as we have shown, will sharpen the 
cosmological limit to the number of light particle species. 
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Figure 1: Big-bang production of the light elements for N v = 2 (short-dashed line), N u = 3 
(solid line), and N u = 4 (long-dashed line). Note that for low baryon density the yields of 
both 7 Li and 4 He increase with N v so that 7 Li and 4 He can be used together to limit N v . 
(For clarity, the theoretical uncertainty in tne light-element production is not shown; see 
e.g., Ref. I.) 




Figure 2: Marginal likelihood for N = N v — AY/0. 016 with different Bayesian priors for 
the primeval deuterium abundance: (D/H)p < 1.0 (solid line); [(D + 3 He)/H]p < 2 x 10~ 4 
(short-dashed line); extreme model of 3 He chemical evolution (from Ref. |22|| ) (long-dashed 
line); (D/H)p = (2.5±0.5) x 10~ 5 (dashed-dotted line). In each case we have assumed the 7 Li 
abundance that results in the least stringent limit to N. The fact that N = 3 is well within 
the 95% credibility interval is indicative of the consistency of big-bang nucleosynthesis with 
three massless neutrino species. 
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